Optical far-infrared properties of graphene monolayer and multilayers 
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We analyze the features of the graphene mono- and multilayer reflectance in the far-infrared region 
as a function of frequency, temperature, and carrier density taking the intraband conductance and 
the interband electron absorbtion into account. The dispersion of plasmon mode of the multilayers is 
calculated using Maxwell's equations with the influence of retardation included. At low temperatures 
and high electron densities, the reflectance of multilayers as a function of frequency has the sharp 
downfall and the subsequent deep well due to the threshold of electron interband absorbtion. 
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Monolayer and bilayer graphenesii^ii^ are gapless two- 
dimensional (2D) semiconductora^i^ whereas its 3D 
predecessor — graphite is a semimetalii^ii^. Hence the 
dimensionality effects for the unique substance can be 
studiedly. Monolayer graphene has a very simple elec- 
tron band structure. Near the energy e = 0, the energy 
bands are cones ei,2(p) = ^vp at the K points in the 
2D Brillouin zone with the constant velocity parameter 
V = 10^ cm/s. Such a degeneration is conditioned by 
symmetry because the small group Csv of the K points 
has two-dimensional representation. 

While the carrier concentration is decreasing in the 
field gate experiment, the graphene conductivity at low 
temperatures goes to the finite minimal valuesi>^. Much 
theoretical effort o^^'^^i^^'^'* have been devoted to eval- 
uate the minimal conductivity in different approaches. 
Theoreticali^i^i^ and experimental researches show that 
the main mechanism of the carrier relaxation is pro- 
vided by the charged impurities and gives the collision 
rate ~ 2iT'^e^nijnp/he'^g£^ where Cg is the dielectric 
constant of graphene, e is the characteristic electron en- 
ergy (of the order of the Fermi energy or temperature), 
and riimp is the density of charged impurities per the 
unit surface. Plasmons in graphene are considered in 
Refs."'^^'^^. The optical visibility of both monolayer and 
bilayer graphene is studied in Ref.^*^ focusing on the role 
of the underlying substrate. 

In the present paper, we analyze the spectroscopy of 
the graphene monolayer and multilayers in the infra-red 
region. In order to calculate the reflection coefficient for 
the multilayers, we follow the method used in Ref.— and 
determine the spectrum of electromagnetic excitations — 
plasmons. We use the appropriate boundary conditions 
at interfaces and the complex conductivity ct as a function 
of frequency w, temperature T, and chemical potential 
11. The chemical potential of ideal pure graphene equals 
to zero at any temperature. With the help of the gate 
voltage, one can control the density and type (n or p) of 
carriers varying their chemical potential. 

The general expression for the conductivity used here 
is obtained in our previous paper— and is valid under a 
restriction that the collision rate of carriers is less than 



the frequency and spatial dispersion of the electric ac 
field, <^ uj, kv. In limiting cases, our result coincides 
with the formulas of JieiM'^^ . For high frequencies, when 
one can also ignore the spacial dispersion of the ac field, 
Lo ^ kv,T~^, the complex conductivity [see Eq. (8) in 
RefiiS] is given by 
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Here, the first term corresponds to the intraband 
electron-photon scattering processes. One can obtain 
it from the Drude-Boltzmann expression (for a case 
1/r = 0) and write explicitly: 
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The second term in Eq. (1), where (5 — > is the in- 
finitesimal quantity determining the bypass around the 
integrand pole, owes its origin to the direct interband 
electron transitions. The real part of this contribution 
is reduced to the expression for the absorbed energy due 
to the interband transitions. Since there is no gap be- 
tween the conduction band and valence band, these two 
terms can compete and the interband contribution be- 
comes larger at high frequencies uj > T, ji. In the op- 
posite case, the intraband contribution plays the leading 
role. 

The difference of the Fermi functions in the second 
integrand equals to 

^ sinh(g/r) 
^ ' cosh{n/T) + cosh(e/T) ' 

Extracting the principal value of the integral, we arrive at 
the integral without singularities and write the interband 
conductivity in the form available for numerical calcula- 
tions: 



+ 00 

m J uj^ ~ ie^ 





(3) 



2 



Here, the first term is given asymptotically by 
G(./2) = ( ^-^(-/^n M<<T, 

\ e(uj - 2^), ^i:$>T, 



(4) 



where the step-function 8{uj — 2fi) expresses the condi- 
tion for the interband electron transitions at low temper- 
atures. The integral in Eq. ([3]) represents the imaginary 
interband correction to the intraband conductivity. 

By using the gate voltage, one can control the den- 
sity of electrons (no) or holes (—no)- Then the chemical 
potential is determined by the condition 



no 



e[/o(e - m) - /o(£ + (5) 



From this expression and Fig. [Ija), one can see that the 
chemical potential goes to zero while the temperature 
increases. 
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FIG. 1: (a) Chemical potential (in K) as a function of tem- 
perature at carrier densities noted at curves, (b) Multilayers 
sample and geometry of wave scattering. 

In order to calculate the graphene reflectance, we apply 
Maxwell's equations 



V(V-E)- V^E = eo^E- 



(6) 



where eo is the ion contribution into the dielectric con- 
stant and j is the conductivity current. We consider the 
case of the p-polarization, when the field E lies in the xz 
plane and the current j has only the in- layer x component 
(see Fig [lb). 

(i) Optics of a monolayer. Consider the graphene 
monolayer at z = with eo = eg deposited on the sub- 
strate (z > 0) with the dielectric constant eo = In 
the vacuum, z < 0, the ac field is given by the sum of 
incident and refiected waves and by the transmitted wave 



in the substrate. In the geometry considered, the current 
in graphene monolayer can be written in the form 



(7) 



Making use of the Fourier transformations with respect 
to the X coordinate, E cx e**^^^, we rewrite the Maxwell 
equations ([6]) as follows 
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FIG. 2: Reflectance from the graphene monolayer with carrier 
densities no — 10* cm~^ (left) and no — 10^° cm~^ (right) 
versus the frequency at temperatures noted at the curves; 
normal incidence. 

The boundary conditions for these equations at z = 
are the continuity of the field component Er^ and the 
jump of the electric-induction z component eEz at sides 
of the monolayer: 

r+O 

fs-Ezl^^+o - ■E^z|z=-o = 47r / p{uj,kx,z)dz . (9) 

The carrier density is connected to the current in Eq. ([7]) 
according to the continuity equation 

p{uj, k^, z) = jx{uj, kx, z)kx/u!. 

Substituting E^ from the second Eq. ([5]) into ©, we find 
the second boundary condition 

e-s dEx 1 dE^ _ 47rcr(w) 

TT^T — U=+o ~ /, \2 ~1 — |2=-o — — ■■ i^x\z=Oi (,J-t)j 

k^ dz \kz) dz tuj 

where 



ks^^esicj/cy^kl kl = ^{cj/c)^~kl 
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Using the boundary conditions, we find the reflection 
amplitude 



1-C 

T — 

1 + C" 

where C = ki[4:Tra{uj)/uj + {es/ks)] 
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FIG. 4: Reflectance from the graphene monolayer with the 
carrier density no — 10^" cm^^ (left) and multilayers with 
no = 10^^ cm~^ in a layer and distance d = 3A between layers 
(right); temperatures are noted at curves, the incidence angle 
is SO''. 



FIG. 3: Reflectance from the graphene monolayer with carrier 
concentrations no = 10^ cm~^ (left) and 10^° cm~^ (right) 
versus temperature for frequencies noted at curves; normal 
incidence. 

The reflection coefficient calculated with the help of 
Eqs. H])-® and pT|) for normal incidence is shown in 
Figs, m and [3] as a function of frequency, temperature, 
and carrier concentration. Notice the different temper- 
ature behavior of reflectance from the samples contain- 
ing the low and high carrier densities. Reflectance for 
the lower carrier density (Fig. [21 left) is larger at higher 
temperatures. This corresponds to the increase in the in- 
traband conductivity (giving the main contribution here) 
with temperature [see Eq. ([2]) for /i < < T] . At the larger 
carrier concentration (Fig. [21 right), the chemical po- 
tential (decreasing with temperature, see Fig. [T^) plays 
the important role appearing in Eq. ([2]) instead of T. 
Therefore, the curve for lowest temperature T = 10 K 
is between the curves for T = 100 K and 300 K in this 
case. The temperature dependencies of reflectance are 
not monotonic as clearly seen in Fig. [3l 

The optical properties of the graphene bilayer can be 
considered in a similar way. Here we do not present 
the corresponding results and investigate another 3D 
example — the graphene multilayers. 

(a) Spectroscopy of graphene multilayers. Let the mul- 
tilayers cross the z axis at points z„ — nd, where d is the 
distance between the layers (see Fig. [TJ)) . Such a system 
can be considered as a model of graphite since the dis- 
tance d = 3A in graphite is larger than the interatomic 
distance in the layer. So we describe the carrier interac- 
tion in the presence of ac electric fleld with the help of 



self-consistent Maxwell's equations (O. For the x com- 
ponent of the fleld , they give 

M \ 
^ + fc2-t-2fc,I?^(5(z-nd) Lb^^O, (12) 

n J 

where 2? = 2i7Ta(Lu)ks/egUj . 

For the infinite number of layers in the stack, the so- 
lutions of Eq. (|12p represent two Bloch states 

ei 2(z) = e±*'=-"''{sinfc^(2 - nd) - e'^'''^'^ x 
sin ks[z - {n + I) d]}, nd < z < {n + l)d (13) 

with the quasi-momentum fc^ determined from the dis- 
persion equation 

cos kzd = cos ksd — V sin ksd . (14) 

The dispersion equation describes the electric field excita- 
tions of the system, i.e., plasmons. The quasi-momentum 
kz can be restricted to the Brillouin half-zone < fc^ < 
Tr/d, if the parameter V is real. In the general case, while 
taking the interband absorbtion into account, we fix the 
choice of the eigen-functions in Eq. by the condition 
Im fcz > so that the solution ei decreases in the positive 
direction z. 

Reflectance from the multilayers occupied the semi- 
space z > can be calculated similarly to the reflectance 
of a monolayer. The electric fleld is given by the de- 
creasing solution ei inside the sample and by the sum of 
incident and reflected waves in the vacuum, 2 < 0, with 
the same value of the component k^ ■ 

Using the boundary condition, Eq. PTI)) . with the di- 
electric constant of graphene Sg instead of Ss , we find the 
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reflected amplitude 

i sin(fcs(i) — Z 
i sin{ksd) + Z ' 

where 

kl is the normal component of wavevector in the vac- 
uum and kz is the quasi-momentum determined by the 
dispersion equation (jl4p at fixed values of w and k^. 

In Fig. m the reflection coefRcient calculated for mul- 
tilayers is shown in comparison with the reflection co- 
efficient of the monolayer. The left panel in this figure 
differs from the right one in Fig. [2] only in the incidence 
angle which is now taken to be 80 ° in order to emphasize 
the multilayer features. The main of them is the sharp 
downfall of reflectance at low temperatures (see the right 



panel in Fig. |4]). This is the threshold effect of the di- 
rect interband transitions at w > 2/x, which is sharp when 
temperature T ^ [see Eq. (j4])]. Just after the downfall, 
the reflectance has the deep well which disappears while 
the temperature increases. Then the downfall becomes 
more smooth. Notice that observations of the absorb- 
tion threshold provide a direct method of carrier density 
characterization of graphene hq = (/i/?iw)^/7r. 

In conclusion, we have developed the detailed micro- 
scopic theory of the graphene mono- and multilayer spec- 
troscopy. We have shown that the nonmonotonic temper- 
ature behavior of reflectance from the monolayer in the 
infra-red region is expected. We have argued that at low 
temperatures and high electron densities, the reflectance 
from multilayers has the sharp downfall with the subse- 
quent deep well. They are caused by the direct interband 
electron transitions. 
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